Search for charged Higgs bosons from top quark decays in 
pp coUisions at y/s = 1.96 TeV 
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We report the results of a search for a charged Higgs boson in the decays of top quarks produced 
in pp collisions at a center-of-mass energy of 1.96 TeV. We use a data sample corresponding to an 
integrated luminosity of 193 pb~^ collected by the upgraded Collider Detector at Fermilab. No 
evidence for charged Higgs production is found, allowing 95% C.L. upper limits to be placed on 
BR(t H^b) for different charged Higgs decay scenarios. In addition, we present in the MSSM 
(m^± ,tan j3) plane the first exclusion regions with radiative and Yukawa coupling corrections. 

PACS numbers: 14.65.Ha, 14.80.Cp, 13.85.Qk 



One of the open questions in the standard model (SM) 
of particle physics involves the mechanism of electroweak 
symmetry breaking (EWSB). Within the SM it is postu- 
lated that a single scalar doublet field breaks the sym- 
metry, resulting in a single observable particle of un- 
known mass called the Higgs boson 1]. To date, the 
SM Higgs boson has not been observed, and extensions 



of the SM have been proposed with different Higgs phe- 
nomenologies. The simplest extension of the SM Higgs 
sector introduces another Higgs doublet. In two-Higgs 
doublet models EWSB results in five Higgs bosons, three 
of which are neutral {h^,H^,A^) and two of which are 
charged (H^). The minimal supersymmetric extension 
of the SM (MSSM) includes a two-Higgs doublet sector 
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in which one doublet couples to the up-type quarks and 
neutrinos, and the other to the down-type quarks and 
charged leptons 0]. The observation of a charged Higgs 
boson would provide unambiguous evidence of a Higgs 
sector richer than that predicted by the SM. 

At the Tevatron, direct production of single charged 
Higgs is expected to be negligible, and the direct pro- 
duction of H^H~ via the weak interaction is expected 
to have a relatively small cross section on the order of 
0.1 pb 0]. The production of ti pairs, with a theoret- 
ical production cross section of 6.7^09 pb |3j for rrit 
=175 GeV/c^, may offer another source of charged Higgs 
production. If kinematically allowed, the top quark can 
decay to H^b, competing with the SM decay t W~^b. 
This mechanism can provide a larger production rate of 
charged Higgs and offers a much cleaner signature than 
that of direct production. 

Previous searches for the charged Higgs boson have 
been performed at y/s = 1.8 TeV in the -|-^T+jets-|-£ 
channels, where the missing energy is defined in 0, 
Th denotes a tau lepton which decays to hadrons, and 
where £ = e or /Lt in Q and £ = e,/i or Th in In 
the framework of the tauonic Higgs model, in which the 
charged Higgs decays exclusively to fi/, these searches set 
limits directly on BR(t H^h) based on the measured 
production rate. These results are then translated into 
limits on tan/3, the ratio of vacuum expectation values of 
the two Higgs doublets. Other searches obtained limits in 
the (m//±,tan/3) plane, assuming that the charged Higgs 
decays only to fv in the -^T+jets-|-T;i channel in j^, and 
assuming that the charged Higgs decays to fv, cs, and 
t*b {-^ W+bb) in the f]T+ jets + i channel with ^ = e or 
/i in [loj . where t* is a virtual top quark. 

These limits utilize tree-level MSSM predictions of the 
t H^b and charged Higgs branching fraction as a 
function of tan p. It is now known that higher-order ra- 
diative corrections significantly modify these predictions. 
The corrections strongly depend on the parameters of the 
model and are particularly large at high values of tan (3 
[l]|. In addition, it is also predicted that in the low tan (3 
region, the charged Higgs has a sizable branching fraction 
to W+h°. 

CDF has recently reported measurements of the ti pro- 
duction cross section in the i -\- fx + jets -I- X channels, 
where £ = e, fj, and where X — £ (the "dilepton" channel), 
X = Th ( "lepton-|-tau" ) , X — one or more tagged jets [l^l 
("lepton-|-jets, > 1 tag"), and X = two or more tagged 
jets ("lepton-|-jets, > 2 tags"). These measurements are 
carried out under the assumption BR(t H~^b) = and 
use data samples corresponding to an integrated lumi- 
nosity of up to 193 pb~^ 0, 13 ^M- Each measure- 
ment agrees with the SM tt cross section expectation 
within its uncertainty, providing no evidence for non-SM 
physics. In this analysis, we consider the possibility of 
t — > and recast the cross section results to set limits 
on charged Higgs production. Depending on the top and 



TABLE L Number of events in each exclusive channel from 
non-tt SM background sources, observed in data, and total 
expected assuming = 6.7 pb and BR{t ~* H'^b) = 0. 



Channel 


Background 


Data 


SM-expected 




events 


events 


events 


dilepton 


2.7±0.7 


13 


10.9 ± 1.4 


lepton-l-jets, = 1 tag 


21.8 ±3.0 


49 


54.0 ±4.3 


lepton-l-jets, > 2 tags 


1.3 ±0.3 


8 


10 ± 1 


lepton + tau 


1.3 ±0.2 


2 


2.3 ±0.3 



Higgs branching ratios, the number of expected events in 
these decay channels can show an excess or deficit with 
respect to SM expectations. 

As published, these measurements allow the catego- 
rization of a single event in multiple channels. In this 
analysis extra requirements are applied to each channel 
in order to force the association of every event to a single 
channel. The ti signal acceptance and non-tt SM back- 
ground contribution to each of these exclusive channels 
are recalculated, and the changes from the original cross 
section analyses are found to be mostly negligible. The 
only exception to this is the "> 1 tag" and "> 2 tags" 
lepton+jets channels, where the latter is a proper subset 
of the former. Removal of this 100% overlap changes the 
"> 1 tag" channel to exactly one tag, "= 1 tag" . The 
results for these new exclusive channels in terms of back- 
ground, number of observed events, and number of SM 
expected events are shown in Table 

We assume that the charged Higgs boson can decay 
only to fv, cs, t*b or W'^h°, leading to five possible decay 
modes for a single top quark: (1) t W~^b, (2) t 
H+ fv, (3) t H+b, H+ cs, (4) t H+b, 
H+ t*b, and (5) t H+b, H+ W+h°, h° bb. 
Charge conjugated decays are implied. Allowing for a 
non-zero BR{t — > H+b), the acceptance of the detector 
for a given ti channel k is 

5 

A = ^ B,- Bj ■ Cij.kiTt, , mH± , m^jo), (1) 

ij = l 

where Bi (Bj) represent the branching fractions of the 
top quark (anti-quark) to decay via mode i (j) as listed 
above, and eij^k is the efficiency to detect a ti event whose 
top quarks decay via modes i and j in channel k. 

The efficiencies £ij,k are obtained from Monte Carlo 
(MC) simulation of ti events generated with different 
masses of the top, H+, and The MC generator 

PYTHIA 16] is modified to include the decay H+ t*b 
and is used for the generation of the ti events. The detec- 
tor simulation and reconstruction algorithms for muons, 
electrons, and jets are identical to those used in the SM 
ti cross section measurements for the four channels. MC 
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efficiencies are scaled for known differences between MC 
simulation of the detector response and that observed in 
data. The dependence of the efficiencies on the width 
of the top quark (Ft) and the width of the charged 
Higgs (r^±) is taken into account using the simulated 
tt events. The systematic uncertainties on e^t for the 
process tt ^ W+bW-b are listed in and do 

not differ much for the other possible decay modes. 
The expected number of events in channel k is 



,back 



(2) 



where cr^™'^ is the tt production cross section and p rep- 
resents a generic model from which the nine quantities 
(five BR's, Tt, r^±, mjj± and m^o), needed to calcu- 
late the acceptance Ak, can be derived. Ck is the inte- 
grated luminosity, and n^"'^'^ is the number of expected 
background events in channel k (shown in Table P). We 
assume the inclusion of the Higgs sector does not modify 
the value of the tt production cross section and set it to 
(jP™'^ = 6.7 ± 0.9 pb. 

For each channel a likelihood is constructed based on 
the Poisson probability to observe events when a 
given model predicts fik events. Since the four channels 
were constructed to be mutually exclusive, the product of 
their likelihoods is taken to form a final likelihood. The 
correlations of the efficiencies, backgrounds, and system- 
atic uncertainties between channels are taken into ac- 
count. The posterior probability distribution of the pa- 
rameter of interest is constructed from the likelihood and 
a prior probability density. The posterior probability is 
integrated to determine the excluded values of the pa- 
rameter. 

In the MSSM the nine quantities needed to calculate 
the acceptance are predicted from a specific set of MSSM 
parameters, including TO^f± and tan/?. We use the com- 
putational package CPsuperH 17] to compute all the 
Higgs masses and branching ratios. This program in- 
cludes QCD, SUSY-QCD, and SUSY-EW radiative cor- 
rections up to the two-loop leading logarithms and ap- 
plies these corrections to the top and bottom Yukawa 
couplings. The top branching ratio to charged Higgs is 
computed with the same level of accuracy from custom 
code developed in collaboration with the authors of ■ 
In the context of the MSSM with to^o < m^± , CP- 
SUPErH predicts that the H+ decay to W^+A" is non- 
negligible for masses of below 100 GeV/c^. In this 
case CPsuperH also predicts the mass of the to be 
similar to that of the /i", and we assume the kinematics 
of the decay to W+A° to be identical to that of W+h° 
when the /i° and A^ masses are equal. Thus, we as- 
sign to the decay W^h*^ a branching ratio of 
BR{H+ ^W+h'^) + BR{H+ ^W+A'^), effectively con- 
sidering both decays. 

As an example of how a charged Higgs alters the bal- 
ance between the top decay channels, Fig.^a) shows the 
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FIG. 1: Predictions for m^^ = 120 GeV/c^ and mt = 175 
GeV/c^ as a function of tan /3 for 193 pb"^ The MSSM pa- 
rameters are defined in [T^ and are set to Msvsy = 1000 
GeV/c^ /i = -500 GeV/c^ At = ^ 2000 GeV/c^, 
Ar = 500 GeV/c^ M2 = = Mq = Mu = Md = Me = 
Ml = MsuSY, and Mi = 0.498M2. (a) Expected number of 
events in each of the channels, (b) Posterior probability den- 
sity obtained when the number of observed events is equal to 
that expected from the SM. 



expected number of events in each of the exclusive chan- 
nels as a function of tan/5 for m^fi = 120 GeV/c^. The 
other relevant MSSM parameters are detailed in the cap- 
tion. The figure demonstrates the excess expected in the 
lepton-|-tau channel and the deficit expected in the other 
channels for large tan /3 values. For values of tan /? around 
7 the BR(i goes to zero and the SM expectation 

for the different channels is recovered. The relationship 
between the channels changes with charged Higgs mass. 
Values of tan (i which result in a non-self-consistent Higgs 
sector are reported by CPsuperH and are shown as the 
theoretically inaccessible regions in Fig. ^ 

Figure db) shows the posterior probability obtained 
for the four channels when the number of observed events 
is equal to that expected from the SM. The posterior is 
obtained by means of a flat prior in logio(tan/3). This 
prior allows for a smooth variation of the top and charged 
Higgs branching ratios as a function of log]^Q(tan/3). The 
probability is integrated over its maximum density region 
to obtain expected upper and lower limits on tan /3 at the 
95% confidence level (C.L.) 
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FIG. 2: The MSSM results obtained with 193 ph'^ at CDF. 
The SM-expected exclusion limits are indicated by black solid 
lines and the zblcr confidence band around it is obtained by 
generating pseudo-experiments. The darkest solid region rep- 
resents the area excluded at 95% C.L. The solid lower region is 
the LEP combined results from direct searches [l9| . Other rel- 
evant MSSM parameters are detailed in the caption of Fig. Q 



and 1 and null elsewhere. The 95% C.L. is obtained 
by integrating the posterior over the maximum density 
region. This procedure is repeated for different charged 
Higgs masses. In the region 80 GeV/c^ < TO//± < 160 
GeV/c2 we exclude BR{t H+b) > 0.4 at 95% C.L. 

Finally, in order to reduce the model dependence, we 
place limits on BK{t H^b) that hold for any combi- 
nation of charged Higgs branching ratios. For a specific 
charged Higgs mass we divide each of the five charged 
Higgs branching ratios into 21 bins. This results in 
1771 possible branching ratio combinations subject to 
the relation Xi) = 1. For each combi- 

nation, we obtain a limit on BR(i —^ H'^b) assuming 
BR(/iO ^ bh) = 0.9 and mno = 70 GeV/c^ The least 
restrictive limit is quoted and the analysis is repeated for 
different charged Higgs masses. The results are shown 
in Fig. 01 For BR{t H+b) > 0.9 (hatched region) the 
width of the top is larger than 14 GeV/c^ and the analyt- 
ical corrections to the efficiencies start losing accuracy. 



Using the number of events observed in the data and 
repeating this procedure for different Higgs masses re- 
sults in the exclusion region shown in Fig. |21 We de- 
termine this exclusion region for several sets of bench- 
mark parameters, including the maximal and minimal 
light Higgs mass scenarios described in p'oj . The com- 
plete characterization of these scenarios and their re- 
sults are described in jij]. In all the benchmarks used, 
the low tan (3 region is excluded in a similar region as 
shown in Fig. |21 The high tan /3 exclusion region, how- 
ever, can be significantly reduced and even vanishes, due 
to parameters of particular benchmarks that suppress 
BR(t H+b). The obtained exclusion limits strongly 
depend on the prior probability used. Using a fiat prior 
in tan /3, which is characterized by sudden changes in 
the top and charged Higgs branching ratios, yields sig- 
nificantly different exclusion regions. It is important to 
note that even if all the corrections were turned off, and 
tree-level calculations were used, the results would be sig- 
nificantly stronger than those obtained in (lOl] under the 
same conditions. 

In the high tan/3 region the decay H+ — > fv is ex- 
pected to dominate in a large fraction of the MSSM pa- 
rameter space. In this region the tauonic Higgs model is a 
good approximation, and we explicitly set BR(i?~'" — > fiy) 
— 1 and evaluate the posterior probability as a function 
of BR(t H+b). The value of r^± has little effect on 
the results as width corrections to the efficiency are small; 
we set r^± =1.4 GeV/c^. The width of the top is set to 
Ff = i-B\(t'^g+b) ' value of m/jo is irrelevant in 

this model. We perform the scan in BR(t H+b) from 
to 1. A posterior probability density of BR(i H+b) 
is obtained using a fiat prior that is constant between 



In summary, we have performed a search for a charged 
Higgs boson in top quark decays using measurements of 
the top pair production cross section in four different 
final states and we find no evidence of signal in the re- 
gion 80 GeV/c2 < mH± < 160 GeV/c^ In the con- 
text of the MSSM with full radiative corrections we ex- 
clude the low tan/3 region for all benchmarks in [2ll |. 
The high tan /? region cannot be excluded independent 
of the MSSM parameters. In the tauonic Higgs model, 
in which the charged Higgs decays exclusively to fv, the 
BR(t H+b) is constrained to be less than 0.4 at 95% 
C.L. If no assumption is made on the charged Higgs de- 
cay, the BR(i H+b) is constrained to be less than 0.91 
at 95% C.L. 
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FIG. 3: Results for the charged Higgs branching ratio inde- 
pendent analysis with nit = 175 GeV/c^. The dark solid 
region represents the CDF Run II excluded region in the 
{■mjj±,BK{t — > H'^b)) plane. The expected exclusion lim- 
its are indicated by a black solid line and the ±lcr confidence 
band around it is obtained by generating pseudo-experiments. 
The hatched region of BR{t — > H'^h) > 0.9 indicates that the 
width of the top is larger than 14 GeV/c^. 
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